The probability distribution of a trapped Brownian particle in plane shear flows 
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We investigate the statistical properties of an over-damped Brownian particle that is trapped by a 
harmonic potential and simultaneously exposed to a linear shear flow or to a plane Poiseuille flow. Its 
probability distribution is determined via the corresponding Smoluchowski equation, which is solved 
analytically for a linear shear flow. In the case of a plane Poiseuille flow, analytical approximations 
for the distribution are obtained by a perturbation analysis and are substantiated by numerical 
results. There is a good agreement between the two approaches for a wide range of parameters. 

PACS numbers: 05.10.Gg, 05.40.-a, 83.50.Ax 



The dynamics of Brownian particles in fluids is of cen- 
tral importance in many areas of science 0, Q . There is a 
profound understanding of Brownian motion in quiescent 
fluids, but the situation is different for particles in flows. 
Several statistical properties of free Brownian particles in 
open flows have been investigated, for instance, in terms 
of the corresponding Langcvin and Fokkcr-Planck equa- 
tions Recent works also consider inertia effects on 
diffusion in shear flows (3, HJ and the relation between 
the Gaussian nature of noise and time-reversibility in 
driven systems @ . It has also been shown that shear flow 
causes, compared to a quiescent fluid, additional correla- 
tions of the particle's velocity and positional fluctuations 
[USB- However, the detection of statistical properties 
of free particles in a flow is an intricate issue. 

To overcome part of these problems, three recent works 
focused on the dynamics of Brownian particles trapped 
by harmonic potentials and exposed to shear flows 0- 
UA\ . It was shown that the shear- induced correlations 
between the positional fluctuations of a captured particle 
along orthogonal directions are essentially the same as in 
the free-particle case (lfj| . Furthermore, a surprisingly 
good agreement was found between the predictions and 
the measurements of these correlations [9| . 

In Ref. [l(| the probability distribution of a trapped 
Brownian particle in shear flows was obtained by simu- 
lations of the Langevin equation. In this brief report the 
related Smoluchowski equation for the probability distri- 
bution is presented and for a linear shear flow an exact 
analytical solution is given. For a plane Poiseuille flow 
we determine approximate analytical solutions, which are 
in good agreement for a wide parameter range with nu- 
merical solutions of the Smoluchowski equation and with 
simulations of the Langcvin equation as well. 

We consider a Brownian particle trapped by a har- 
monic potential with the spring constant k at the origin 
of the Cartesian coordinate frame (x,y, z): 



V(f) 



(1) 



The particle is exposed to a flow along the x-direction 
with a ^-dependent magnitude: 



Since the Brownian motion perpendicular to the shear 
plane is decoupled from the one in the shear plane, we 
consider the quasi two-dimensional problem with r = 
(x,y) and V = (<9 2 ,<9 S ). 

For a plane Poiseuille flow the position of the potential 
minimum can be different from the center of the flow. In 
the shifted coordinate frame, (x, y) = (x, y + yo), where 
?/o describes the y-coordinate of the potential minimum, 
the flow profile is given by u(y) = u p (l — y 2 /l 2 )e x with 
the flow velocity u p at its center and the confining plane 
boundaries at y = ±i. If a particle is trapped at yo =^= 0, 
one obtains with y = yo + y the coefficients a = u p (l — 
y 2 /l 2 ), b = -2u p y /l 2 and c = -u p /l 2 in Eq. ©. This 
work focuses on situations where the particle positions 
are sufficiently far away from the boundaries, so that 
hydrodynamic interactions with the wall, as discussed in 
[3( for instance, can be neglected. 

The particle dynamics are determined by thermal mo- 
tion, potential forces and drag forces caused by the flow. 
The thermal motion is characterized by the diffusion con- 
stant D = fcsT/^, which is given by the Einstein relation 
in terms of the temperature T, the Boltzmann constant 
ks and the Stokes friction coefficient £ = 6nr)R, where r\ 
is the viscosity of the fluid and R is the effective radius 
of the particle. The external flow u is the origin of the 
Stokes drag force £u on the point-like particle, which is 
balanced by the restoring force — VV = —kr. 

The diffusion and the two deterministic forces drive the 
probability current i(x,y,t) of the Smoluchowski equa- 
tion (SE) for the particle's positional distribution func- 
tion V(x,y,i) 0: 



diP = -Vj , 

j = -DVV + (u - ^VV)V . 



(3) 
(4) 



With the expressions in Eqs. (JTJ - Q the SE of a 
Brownian particle trapped in a harmonic potential and 
exposed to a shear flow reads: 



d{P =DV 2 V + ^(2 + xd s + ydy) V 
- (a + by + cf)d s P. 



u(r) = (a + by + cy 2 )e x 



(2) 



(5) 

The two spatial coordinates may be rescaled by the 
length S = ^JksT /k, x = Sx and y = 5y, alike the time, 
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FIG. 1. An elliptical contour line of the distribution of a 
trapped particle in a linear shear flow is shown, cf. Eq. (jHJ 
with a = 0, P = 1, as well as the vector field j(x,y). The 
broad arrows indicate the two principal axes of the ellipse. 



t = Q/kt. This results in the dimensionless SE 

d t V = [V 2 + 2 + yd y + (x-a-/3y- jy 2 )d x ] V , (6) 

with the parameters a = aS/D, j3 = bS 2 /D and 7 = 
cS 3 /D describing the flow profile and y G [(— I— yo)/S 7 (I— 
yo)/S\. (3 is the so-called Weissenberg number. We note 
here, that a modified Smoluchowski equation including 
inertia in shear flows is presented in [J, Q . In comparison 
to these works, the presence of the potential (JXJ) ensures 
a stationary solution V(x, y). 

For a uniform flow, i.e., a ^ and /3 = 7 = 0, the static 
solution of Eq. <j6j) is the shifted Boltzmann distribution 



V n {x,y) = P a e' 



1 2 1 



(7) 



where Pq = (2ne a / 2 ) _1 is determined by 
J J dxdyP n (x,y) = 1. A superposition of a uni- 
form flow and a linear shear flow, i.e., a 7^ 0, /3 ^ and 
7 = 0, leads to the expected Gaussian distribution 



Pab(x,y) = P a e 
with the coefficients 
2 
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and the norm P = U yj(3 2 +Ae 2a +4 )J . In the 

case of a finite Weissenberg number (5 the contour lines 
of the probability distribution P a b{x,y) are elliptical, as 
shown in Fig. [T] The probability current j(x,y), indi- 
cated by the vector field in the same figure, characterizes 
the mean particle motion. 

The parameter a, describing the contribution of the 
uniform flow, causes essentially a shift of the distribution 



T-'ab{x,y) in the x-direction. In the following we consider 
the case a = 0, where et 4 = 05 = and the resulting 
distribution function is denoted by V a (x,y). Its ellipti- 
cal contour lines can be characterized by two principal 
axes. The ratio between their lengths as well as the an- 
gle between the major principal axis and the flow direc- 
tion are a function of the Weissenberg number /?. This 
was already discussed in Ref. [loj . where shear- induced 
corrections to the autocorrelations as well as a cross- 
correlation between orthogonal particle fluctuations in 
the shear plane were found. These correlations can also 
be calculated in terms of the probability distribution 
V a {x, y) via the expression (r^rj) = J J dxdyV a {x,y)rirj: 
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(xy) = (yx) = £ (yy) = l. (10) 



After rescaling to dimensional units, x — > x, the results 
given in Ref. [l(| are recovered. 

In a plane Poiseuille flow all three parameters a, (3 and 
7 in the SE © may be non-zero and no exact analyti- 
cal solution was found in this case. Similar to Q, the 
probability distribution is calculated pcrturbatively and 
compared with numerical solutions of Eq. ([6]) . 

First we consider a Brownian particle which is trapped 
at the center of a plane Poiseuille flow, with yo — 0. In 
this case of vanishing /3, the parameter a = u p 8^/(kBT) 
describes the ratio between the drag force imposed by 
the flow and the potential force on the particle. If a is 
fixed, the parameter 7 = ~a(S/l) 2 depends on the ratio 
between the two characteristic lengths / and 6. Note, the 
hydrodynamic interactions between the particle and the 
walls arc small in the range 6/1 < 1/2. 

Our ansatz for the perturbation expansion of the solu- 
tion of Eq. <j6j) up to the second order in 7 reads 



V(x,y) =V n (x,y) e T/i(*,i0+7V(*.v) 



(11) 



with the two polynomials 

fi{x,y) = hx + b 2 y 2 + b 3 xy 2 , (12a) 
gi(x,y) = (c 1 x + c 2 x 2 )(l + y 2 ) + c 3 y 2 + c 4 y 4 . (12b) 

The SE ([6]) may then be rewritten: 

{Po{x,y) + jpi(x,y) + j 2 P2(x,y)) V(x,y) = 0. (13) 

Since 7 is an arbitrary, but small number, the polyno- 
mials po, 1,2(2, y) in Eq. (|13[) have to vanish separately 
According to Eq. @ the condition Po(x,y) — is au- 
tomatically fulfilled. The second condition, pi(x,y) = 0, 
provides the first-order coefficients 
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whereas the third condition, p 2 {x, y) 
coefficients at C(7 2 ): 
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(14) 

0, determines the 
04 = ^. (15) 
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With the potential minimum off the center of the 
Poiseuille flow, one has yo 7^ 0, a finite Weissenberg num- 
ber P ^ and no longer ±y symmetry. Again we use an 
ansatz of the form: 



V{x,y) = T ab {x 1 y)e~ <h(x > v)+ ^ 92(x > v) 



(16) 



Due to the loss of the ±y symmetry, the polynomial 
f2(x,y) has nine different contributions: 



f2(x, y) =dix + d 2 y + d 3 (x - 3a)x 
+ d$(x - 2a)xy + d e y 3 . 



d/±(x — a)y 2 



(17) 



The condition that the 0(7)-terms in Eq. ([6]) have to van- 
ish determines the coefficients in the ansatz p7|) . With 
B = [3 2 + 4 and E = 2/(9B 3 ) they are given by 



di = 4£ ((12 - /3 2 )B + 16a 2 /3 2 ) , 
d 2 = -2E/3 ((20 + f3 2 )B + 8a 2 {/3 2 - 4)) , 
64 

d 3 = —E/3 2 , d 4 = 2£(3/3 4 - 8/3 2 + 48), 

d 5 = 16£/3(4-/3 2 ), 

-Ef3 



-(5/3 4 + 24/3 2 + 144). 



(18) 



In the limit /3 — >■ the coefficients ([14]) are recovered. The 
polynomial g2(x,y) for the order j 2 includes 14 lengthy 
contributions, which we do not list here. 

In order to estimate the validity range of the approx- 
imations presented above, we determine stationary solu- 
tions of the SE. ([6]) numerically. To this end a simple 
Jacobi-relaxation method, cf. Ref. [l3j], or a direct inte- 
gration of the rescaled Eq. ^ is sufficient. The proba- 
bility current is determined via: 



-VP 







(19) 



Fig. [2] shows the contour lines of the numerically ob- 
tained probability distribution V{x,y) in the case with 
the minimum of the trapping potential at the center of 
a plane Poiseuille flow. This distribution has similarities 
with the parachute or bullet-like shape of vesicles in a 
Poiseuille flow [I3 - [l6j . In comparison to similar numer- 
ical results, obtained via the related Langevin equation 
and presented in Ref. fioj , we additionally show the prob- 
ability current j(x, y). As indicated in Fig. [2J the vector 
field i{x,y) includes two counter-rotating vortices that 
are symmetric with respect to the x-axis. 

In Fig. [3] contour lines of the numerically obtained 
probability distribution are compared with those ob- 
tained from the perturbative solution (jlip for the pa- 
rameters a = 6 and 7 = — a(8/l) 2 = —3/8. In spite of 
this rather large |7|-value, beyond the expected validity 
range of the perturbation expansion, the differences be- 
tween both solutions in Fig. [3] are surprisingly small. As 
expected, for decreasing values of I7I these differences be- 
come even smaller, but the analytical formula pip may 
be useful for fitting experimental data up to I7I ~ 3/8. 




FIG. 2. Contour lines of the numerical solution of Eq. ([6]) as 
well as the vector field j(x,y) are shown for the case where 
the potential minimum is at the center of a plane Poiseuille 
flow. Parameters: a — 8, 7 = — a(5/l) 2 — —1/2. 
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FIG. 3. Comparison of the numerically determined distribu- 
tion (solid lines) with the analytical approximations, given by 
Eq. (dashed lines) . The contour lines of both types are 
plotted at different heights in order to distinguish them to 
compare their shapes. Parameters: a = 6, 7 = —3/8. 



The symmetrical shape of the distribution is lost if the 
particle is trapped away from the center of the plane 
Poiseuille flow. With increasing values of yo and /?, 
the shape of the probability distribution deforms from 
a parachute or bullet towards an ellipse. Simultaneously 
one vortex of the probability current is enhanced while 
the other one is weakened. For the values yo = 1/4, 
5/1 = 1/4 and a = 7, with (3 = -2a/15 and 7 = -a/15, 
the contour lines of V(x, y) as well as the current j(x,y) 
are displayed in Fig. 3] In this case the lower vortex in 
j(x,y) has already vanished and the distribution V{x,y) 
shares similarities with the slipper shape of vesicles in 
capillary flows [T^j- Again the results are in good agree- 
ment with the simulations in Ref. [To| . 

In conclusion, we investigated the positional distribu- 
tion of a Brownian particle, which is trapped by a har- 
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FIG. 4. Contour lines of V(x,y) and vector field ](x,y) for 
the case where the potential minimum is off the center of the 
Poiseuille flow. Parameters: a — 7, /3 — —14/15,7 — —7/15. 



monic potential and simultaneously exposed to different 
shear flows. A complete analytical solution of the cor- 
responding Smoluchowski equation is given in the case 
of a linear shear flow. For a plane Poiseuille flow, we 
presented approximate analytical formulas, which are in 
good agreement with numerical solutions for a wide range 
of parameters. Some of our results confirm earlier ones 
obtained in Ref. on the basis of simulations of the 
related Langevin equation. 

Our predictions of the particle's probability distribu- 
tion in a Poiseuille flow may be measured in an experi- 



ment similar to that in Ref. In this work, a micron- 
sized polystyrene bead was trapped by an optical tweezer 
and the fluctuating particle positions were recorded by a 
high speed camera in a stroboscopic manner. For the 
positions of a trapped particle in Poiseuille flow one ex- 
pects probability distributions as predicted in this work. 
The presented analytical approximations for the particle 
distribution may be useful for fitting the experimental 
data. 

Experiments with particles placed near the center of 
a Poiseuille flow, where the flow velocity takes its max- 
imum value, may require large laser intensities in order 
to keep the particles trapped by the potential. This con- 
straint reduces the flexibility for variations of the typi- 
cal length scale S of the particle's positional fluctuations. 
Since the flow profile determines the shape of the particle 
distribution and not the flow velocity near the potential 
minimum, one may reduce the drag force by moving the 
trap along with the flow. 

The statistical properties of trapped Brownian parti- 
cles, as investigated in 0, [l(| EH], share similarities with 
those of tethered polymers exposed to uniform flows flfjl - 
|22| or to shear flows (2^|. The related theoretical explo- 
rations are mainly based on Brownian dynamics simula- 
tions. An interesting question is, how far can the analysis 
described here be applied to tethered bead-spring models 
in shear flows? Do such investigations exhibit temporal 
oscillations as found for instance for deterministic models 

H? 
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